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Abstract. The propagation of plane waves in a Faraday chiral medium is investigated. Condi- 
tions for the phase velocity to be directed opposite to the direction of power flow are derived for 
propagation in an arbitrary direction; simplified conditions which apply to propagation parallel to 
the distinguished axis are also established. These negative phase-velocity conditions are explored 
numerically using a representative Faraday chiral medium, arising from the homogenization of an 
isotropic chiral medium and a magnetically biased ferrite. It is demonstrated that the phase velocity 
may be directed opposite to power flow, provided that the gyrotropic parameter of the ferrite com- 
ponent medium is sufficiently large compared with the corresponding nongyrotropic permeability 
parameters. 

PACS number(s): 41.20.Jb, 42.25.Bs, 83.80.Ab 



1 Introduction 

Homogeneous mediums which support the propagation of waves with phase velocity directed op- 
posite to the direction of power flow have attracted much attention lately [1, 2, 3]. The archetypal 
example of such a medium is the isotropic dielectric-magnetic medium with simultaneously nega- 
tive real permittivity and negative real permeability scalars, as first described Veselago in the late 
1960s [4]. A range of exotic and potentially useful electromagnetic phenomenons, such as negative 
refraction, inverse Doppler shift, and inverse Cerenkov radiation, were predicted for this type of 
medium [3, 4]. Recent experimental studies involving the microwave illumination of certain com- 
posite metamaterials [5, 6] — which followed on from earlier works of Pendry et al. [7, 8] — are 
supportive of Veselago's predictions and have prompted an intensification of interest in this area. 
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In particular, a general condition — applicable to dissipative isotropic dielectric-magnetic medi- 
ums — has been derived for the phase velocity to be directed opposite to power flow [9]; and this 
condition shows that the real parts of both the permittivity and the permeability scalars do not 
have to be negative. 

A consensus has yet to be reached on terminology. For the present purposes, a medium support- 
ing wave propagation with phase velocity directed opposite to power flow is most aptly referred to 
as a negative phase-velocity medium. However, the reader is alerted that alternative terms, such as 
left-handed material [3], backward medium [10], double-negative medium [1], and negative-index 
medium [11], are also in circulation. A discussion of this issue is available elsewhere [12]. 

The scope for the phase velocity to be directed opposite to power flow may be greatly ex- 
tended by considering non-isotropic mediums, as has been indicated by considerations of uniaxial 
dielectric-magnetic mediums [10, 13, 14]. The focus of the present communication is the propa- 
gation of negative phase-velocity plane waves in Faraday chiral mediums (FCMs) [15, 16]. These 
mediums combine natural optical activity — as exhibited by isotropic chiral mediums [17] — with 
Faraday rotation — as exhibited by gyrotropic mediums [18, 19, 20]. A FCM may be theoretically 
conceptualized as a homogenized composite medium (HCM) arising from the blending together of 
an isotropic chiral medium with either a magnetically biased ferrite [21] or a magnetically biased 
plasma [22]. The HCM component mediums are envisioned as random particulate distributions. 
The homogenization process is justified provided that the particulate length scales in the mixture 
of components are small compared with electromagnetic wavelengths. A vast literature on the esti- 
mation of constitutive parameters of HCMs exists; see Refs. [23, 24], for example. The constitutive 
relations of FCMs have been rigorously established for some time [16], although inappropriate use 
still occurs [25]. 

In the following sections, wavenumbers and corresponding electric field phasors are delineated 
from eigenvalue/vector analysis for planewave propagation in an arbitrary direction. Simpli- 
fied expressions for these quantities are derived for propagation parallel to the biasing (quasi)- 
magnetostatic field [19, p. 71]. A general condition for the phase velocity to be directed opposite 
to power flow is established. The theoretical analysis is illustrated by means of a representative 
numerical example: the constitutive parameters of FCMs arising from a specific homogenization 
scenario are estimated and then used to explore the wave propagation characteristics. 

As regards notation, vectors are underlined whereas dyadics are double underlined. All electro- 
magnetic field phasors and constitutive parameters depend implicitly on the circular frequency 
ijj of the electromagnetic field. Unit vectors are denoted by the superscript " symbol, while 
I_ = xx + yy + z_z_ is the identity dyadic. The complex conjugate of a quantity q is written as 
q*; the real part of q is written as Re{g}. The free-space (i.e., vacuum) wavenumber is denoted 
by k = io^e fj, where e and fi are the permittivity and permeability of free space, respectively; 
and r] = yTW^o represents the intrinsic impedance of free space. 
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2 Analysis 



2.1 Preliminaries 

The propagation of plane waves with field phasors 

= Ko exp(ik ku • r) 
H(r) = H_ exp(ik ku • r) 



(1) 



in a FCM is considered. Such a medium is characterized by the frequency-domain constitutive 
relations [16] 

D(r) = g • E(r) + £ • H(r) \ 
B{r) = -£ - E{r) + fi • H(r) ] ' 



(2) 



with constitutive dyadics 

e = e [ e / — ie g z_ x / + ( e 2 — e ) ||] 

| = z^e^ [CZ-^lxI+(^-0 11] (3) 

Thus, the distinguished axis of the FCM is chosen to be the z axis. For FCMs which arise as 
HCMs, the gyrotropic parameters e g , £ g and jig in (3) develop due to the gyrotropic properties of 
the ferrite or plasma component mediums. Parenthetically, it is remarked that more general FCMs 
can develop through the homogenization of component mediums based on nonspherical particulate 
geometries [22, 26]. 

In general, the relative wavenumber k in (1) is complex valued; i.e., 

k = ~k R + ik I , fefe/el). (4) 
It may be calculated from the planewave dispersion relation 

det L(ik ku) =0, (5) 
which arises from the vector Helmholtz equation 

UV)-E(r) = 0, (6) 

wherein 

UV) = (V x l + iuj^j - ff 1 - (V x l + iio^j -uj 2 g. (7) 

Of particular interest is the orientation of the phase velocity, as specified by the direction of 
kjiu, relative to the direction of power flow given by the time-averaged Poynting vector P(r) = 
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iRe [E_(r) x H_*(r)]. The combination of the constitutive relations (2) with the source-free Maxwell 
curl postulates 

V x E(r) = icoB(r) 

V x H(r) = -iuD(r) 

yields 

P(r) = - exp (^—2k kiu • rj Re ji? x 
for plane waves (1). 

In the remainder of this section, the quantity kn u • P(r) is derived for planewave propagation 
in an arbitrary direction; without loss of generality, u is taken to lie in the xz plane (i.e., u = 
xsin9 + z_ cosO). Further manipulations reveal the simple form Uru • P(r) adopts for propagation 
along the FCM distinguished axis (i.e., u = z). 




Qf 1 )* ' (Vw k*uxEZ+i* *^o)]} 



(9) 



2.2 Propagation in the xz plane 

For u = x sin9 + z cos 9, the dispersion relation (5) may be represented by the quartic polynomial 



with coefficients 

CL4 = 
(13 = 

a 2 = 



ai 
a 



a 4 k 4 + a 3 k 3 + a 2 k 2 + ajt + a = , (10) 

= (e sin 2 9 + e z cos 2 9) (fi sin 2 9 + n z cos 2 9) - (£ sin 2 9 + £ z cos 2 9) 2 , (11) 
= 2 cos 9{ sin 2 9 [n g (e£ 2 - e z £) + e g (/x£ 2 - /j, z £) + £ g (fie z + efi z - 2££ g )] 

+2cos 2 ^ 9 (e^-£ 2 )}, (12) 

= sin 2 e[nn z (e 2 - e 2 ) + (f + £ 2 ) (//e* + efi z ) - 2£ [6 (? 2 - £ 2 ) + H g e z S g ] 
-2e g [£ g (n g $, - ^g) + nzttg] - e [e z (fi 2 - fi 2 g ) + 2£ 2 (//£ - fi g Q] } 
+2 cos 2 9 (e zt i z - f z ) (3£ 2 - f - e gf i g - e/i) , (13) 
= Acos9{e z n z - £ 2 Z ) £ (e 9 ^ + e/x s ) + £ 9 (£ 2 - £ 2 - e/i - €g fi g ) , (14) 

= {e zllz - e z ) [ {e 2 ~ e 2 g ) (^ 2 - ^) + {f g - f f ~ 2 {f g + f) (^ + e g ^ g ) 

+^g(^g + ^ g )]. (15) 



Hence, four relative wavenumbers k = Ki, k«, Km and Ki v may be extracted — either algebraically 
or numerically [27] — as the roots of (10). 

Upon substituting u = x sin 9 + z cos 9 into (9) and combining with (3), the component of P(r) 
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aligned with u emerges straightforwardly as 
u ' P(z) = 7T - exp (-2k k!U • r) Re J -\ (k* sine\E 0y \ 2 - iC E OyKz) sin # 



Mi 



fc* [M* {\E 0x \ 2 + l^o/) + ^ {E 0x E* y - E 0y E* x )] cos 2 



+^*\E 0z \ 2 sin 2 - [/x* (E 0z E* x + So^oJ + ^ {E 0z E* y - E 0y E* z )] sm6 cos 
+ " ^t) [{\E 0x \ 2 + \E 0y \ 2 ) cosO- E 0z E* 0x sm9] 

-i (/x*r - M^;) [(^oa^oj, - E* 0x E 0y ) cos9- E 0z E* 0y sine] 



wherein (E 0x , E 0y , E 0z ) = E^. 
Let the quantity 

w = 2r/ exp (2k kiu • |£th/r 2 fc/?ii * 
be introduced such that the fulfilment of the negative phase-velocity condition 

k R u'P(r) < 

is signaled by w < 0. 

Substitution of (16) in (17) yields the expression 

= k R Re I \ (k* sin - sin 

1 



(^)2-(^)2 



fc* [/x* (|a| 2 + 1) + ^* (a - a*)] cos 2 



+//|/?| 2 sin 2 6>- [fi* (a* f3 + a(3*) + ifi g (f3 - (3*)] sin^cos^ 
+ (m*^ - V* g C) [(l"| 2 + 1) cosfl - a*/? sin^] 
-i (//f - ^) [(« - «*) cose- (3 sine] 
The ratios of electric field components 

a = Eo x /Eoy | 



(16) 



(17) 
(18) 



(19) 



(20) 



P — E 0z /E y J 

in (19) are derived as follows: As a function of e, the dyadic operator L of (7) has the form 
L = iK{[L] u xx+[L] 22 yy+[L] 33 zz+[L] u (xy-yx) 

+ [L] 13 (xz + zx) + [L] 2Z (yz-zy) }, (21) 



with components 



ik] 22 
ik] 



23 



e + 



2^er-^(e 2 + r 2 ) 



h 2 -14 



~k 2 sm 2 e -n{e + r 2 ) 



h 2 -14 



e, - 



i\ e + 



Hz 

g _ fik 2 sin 2 9 

^ - Hg 



= I 



^—^ — k sin 9 , 
H 2 - fi 2 



\ ll 2 - p 2 



Hz 



k sin 9 , 



where r = £ 9 + A; cos#. It follows from the vector Helmholtz equation (6) that 



a 







IkUlkU + lkUlk 



123 



[^13[^13-^]n[^33 
[k] 12 ik] 23 -[k} 13 [k] 22 



> . 



ikUk^ + ikUk] 



33 > 



(22) 
(23) 
(24) 
(25) 
(26) 
(27) 



(28) 



2.3 Propagation along the z axis 

The results of the preceding analysis simplify considerably for planewave propagation along the z 
axis (i.e., 9 = 0). The quartic dispersion relation (10) yields the four relative wavenumbers 



Ki = V e + <W^ + Vg -£-£g 
Kii = ~y/ e + e gVf i + % 
l*<>iii = y — e g\> 1 H ~ Hg + £ ~~ £3 
*iv = ~V e ~ e gV^ - Vg + £~£g 



(29) 



and (19) reduces to 

(|a| 2 + 1) ( fc* M * - + M*^) + i (a - a*) ( " H*C + 



u> = kji Re 



(/,*) 2 - (^*) 2 



(30) 



Since the dyadic operator components [£] 13 and [£] 23 are null-valued for 9 = 0, the electric field 
ratios are given as 

<* = -ik] 12 /lk]n 
= 



(31) 



Note that a further consequence of [Ll] 13 = [ Ll ] 23 = is that the time-averaged Poynting vector 
is parallel to the z axis. 
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By substituting (29) into (22) and (25), the ratio a emerges as 

a — < * ^ = (32) 

Hence, negative-phase velocity propagation along the z axis occurs provided w < where 



[ i e * -)_ e * I 

w = = 2 Re{ v / e + e 9V / ^ + ^ 3 - £ - £ g } Re < y \ > for fc = 

[ V A 4 + A* 9 J 

(33) 

to = wu = 2 Re { Je + e g -vAt + /x fl + £ + £ fl [ Re { v r —^ = } for k = ku, 



= w u = 2Re {^TT 9 ^JTTjr g + e + £ g } Re 



w = = 2 Re { ^/e- e gy /fj,- \i g + £ - £ g } Re I V = = \ \ for k = Km, 

I V A* A* fl I 



= = 2 Re { y'e - e gy /n - Hg -£ + £<?} Re { , = _ = } for fe = 

I V A* Mflr J 



(34) 



(35) 



(36) 



3 Numerical results 



In order to further examine the negative phase-velocity conditions derived in Sections 2.2 and 2.3, 
let us consider a Faraday chiral medium (FCM) produced by mixing (a) an isotropic chiral medium 
described by the constitutive relations [17] 

D = e e a E + i^r £ a H ) 
B = -i^H ) £ a E + ^ a H J 

and (b) a magnetically biased ferrite described by the constitutive relations [19, Ch. 7] 

D = e e b E 

B = fi [fi b l-ifi b g zxl+(fi b z -fi b )zz] -H 

Both component mediums are envisioned as random distributions of electrically small, spherical 
particles. The resulting homogenized composite medium (HCM) is a FCM characterized by the 
constitutive dyadics 

fCM = ^ [ e HCM ^ _ le HCM 1 x £ + ( ^CM _ ^HCM ) ^ ] 

f CM = i^[ £ HCM L ~ i$ CM z x I + ( £»™ - £ HCM )zz] . (39) 

ffCM = Mq [ ^HCM j _ ^HCM zXl+[ ^ CM - fi HCM ) zz] J 




7 



Incidentally, a FCM with constitutive dyadics of the form (39) may also be developed via the 
homogenization of an isotropic chiral medium and a magnetically biased plasma [22]. 

The constitutive dyadics e HCM , £ HCM and \x HCM are estimated using the Bruggeman homoge- 
nization formalism for a representative example. Comprehensive details of the Bruggeman formal- 
ism [28, 24] and its implementation in the context of FCMs [21, 22, 29] are available elsewhere. 
Initially, we restrict our attention to nondissipative FCMs; the influence of dissipation is considered 
later in this section. 



3.1 Nondissipative FCMs 

The parameter values selected for nondissipative component mediums are as follows: 

e° = 3.2, £ a = 2.4, fi a = 2; e b = 2.2, / = 3.5, & = 1, ^G[0,4]. 

The permeability parameters for component medium b may be viewed in terms of the semi-classical 
ferrite model as 

H° = 1 + — f 

b <pn }, (10) 



u b 

Hz 



wherein to is the Larmor precessional frequency of spinning electrons and uo m is the saturated 
magnetization frequency [18, 19]. Thus, the parameter values fi b = 3.5 and /j, b 6 [0,4] correspond 
to the relative frequency range (uj /uj) £ [0.625, oo). 

Let f a denote the volume fraction of the isotropic chiral component medium a. In figure 1, 
the estimated constitutive parameters of the HCM are plotted as functions of f a for fi b = 4. The 
uniaxial and gyrotropic characteristics of a FCM are clearly reflected by the constituents of the 
permeability dyadic /j, HCM and the magnetoelectric dyadic £ HCM . In contrast, the HCM is close 
to being isotropic with respect to its dielectric properties. Significantly, eight of the nine scalars 
appearing in (39) are positive, while e g M is negative only for f a < 0.32; however, M | << 
\e HCM \ and |ef CM \ « |ef CM \ for all values of f a G [0, 1]. 

The permeability parameters \i HCM and fig are equal at f a w 0.25, it being clear from 
the right side of (19) that this equality has an important bearing on the stability of w. Further 
calculations with fi b = 2 and fi b = 3 have confirmed that fi HCM / l^ CM for all volume fractions 
f a G [0, 1]. This matter is pursued in figure 2 where the estimated constitutive parameters of the 
HCM are graphed as functions of fi b for f a = 0.35. The HCM gyrotropic parameters ^ and 
nf CM are observed to increase steadily as /j, b increases; e^ CM , i^ CM and [i^ CM all vanish in the 
limit /j, b — > 0. Also, as n g increases, the degree of uniaxiality (with respect to the z axis) increases 
for f CM but decreases for fi HCM . 

The relative wavenumbers k = Ki-i v for propagation along the z axis, as specified in (29), are 
displayed in figure 3 as functions of f a , for fi b = 2,3 and 4. The relative wavenumbers K{ > and 
Ku < for all f a € [0, 1] for fi b = 2,3 and 4. Similarly, for fi b = 2 and 3, the relative numbers 
Km > and Ki v < 0. 
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However, the equality [i HCM = \i I g iCM , which occurs at f a « 0.25 for ^ = 4, results in and 
Ki V acquiring nonzero imaginary parts as f a falls below 0.25 for fi b = 4. Only the real parts of these 
complex-valued relative wavenumbers are plotted in figure 3. 

Observe that Kj, Km and Ki V > in figure 3, whereas ku < in the volume fraction range 
0.25 < f a < 0.42 with fi b = 4. Furthermore, Re {Km } and Re } > while Ka < for 

f a < 0.25 with n g = 4. In the limit f a — > 0, the relative wavenumbers — > ±v / ?y^ fe ± (i.e., 

the relative wavenumbers of a ferrite biased along the z axis [19]). Also, as / a — > 1, the relative 
wavenumbers Ki-i v — ► iy^/fi^" (i.e., the relative wavenumbers of an isotropic chiral medium). 

The values of u; corresponding to the relative wavenumbers Ki—i v of figure 3, namely Wi—i v , are 
plotted against f a in figure 4 for [i b = 2,3 and 4. The quantities Wi_m > for all volume fractions 
f a € [0,1] with fi b = 2,3 and 4. Thus, for the relative wavenumbers Ki-m, power flows in the 
same direction as the phase velocity. This is the case regardless of whether the phase velocity is 
directed along the positive z axis (as in modes Kj and Km) or directed along the negative z axis (as 
in mode Ka). Both wm 

and Wi v are null valued in those regions where the corresponding relative 
wavenumbers, Km and K{ v , respectively, have nonzero imaginary parts. In addition, wm — ► °o and 
Wiv - ► — o° i n the vicinity of f a = 0.25 for fi b = 4. 

Significantly, wi v < for fi b = 4 at volume fractions f a G (0.25, 0.42) in figure 4. This means that 
the negative phase-velocity condition then holds in the chosen FCM which has been conceptualized 
as a homogenized composite medium. 

In figure 5, the relative wavenumbers «j and Km for 8 = tt/2 (i.e., propagation along the x axis) 
are plotted against f a for jji h g = 2, 3 and 4. The graphs of Ka and Kj„ need not be presented since 
Ki = —Ka and Km = —Hiv 1 For all / a G [0, 1] with fi b = 2 and 3, the relative wavenumbers «j > 
and Kjjj < 0. Similarly, Kj > for /j, b = 4. However, when [i b = 4, it is found that Km < for 
f a > 0.42 but Km possesses a nonzero imaginary part and Re {Km} = for f a < 0.42. In the limit 

f a — > 0, the relative wavenumbers — > ±-y/ e b / n b \J (ix b ) 2 — (/i^) 2 and ±-y/ e b i4 (i- e -) t ne relative 
wavenumbers of a ferrite biased along the x axis [19]). Also, as f a — > 1, the relative wavenumbers 
— ► ±^e a [i a ± ^ a (i.e., the relative wavenumbers of an isotropic chiral medium). 

Figure 6 shows the plots of corresponding to the relative wavenumbers Ki^u of figure 5. 
The graphs of w = wu^ v corresponding to the relative waveumbers Ku^ v are not displayed since 
the equalities Wi = wa and wm = Wi v hold for 9 = tt/2 — as may be inferred from (19)-(28). The 
quantities w^m > at all volume fractions f a € [0, 1] with fi b = 2, 3 and 4. As remarked earlier for 
Ki-m propagation along the z axis, here we have that power flows in the same direction as the phase 
velocity, regardless of whether the phase velocity is directed along the positive x axis (mode «$) or 
along the negative x axis (mode Km)- Furthermore, it is found that wm = in the region where 
the corresponding relative wavenumber Km is purely imaginary (i.e., for f a < 0.42 with /j, b = 4). 

1 When 6 = ir/2, the dispersion relation (10) reduces to a quadratic polynomial in k 2 . 
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3.2 Dissipative FCMs 



The scope of these numerical investigations is now broadened by considering (i) the effects of 
dissipation or loss; and (ii) propagation in an arbitrary direction. Let a small amount of loss be 
incorporated into component medium b by selecting the constitutive parameters of the component 
mediums as 

e a = 3.2, £ a = 2.4, fi a = 2; e b = 2.2 + i5, n b = 3.5 + iS, /4 = l + i0.5<5, fi b g = 4 + i25, 

where the dissipation parameter 5 G [0,0.2]. We focus attention on the region of negative phase- 
velocity propagation along the z axis with relative wavenumber as illustrated by wi v < at 
0.25 < f a < 0.42 in figure 4. 

Real parts of the relative wavenumber Ki v , calculated at the volume fraction f a = 0.35 with 
5 = 0,0.1 and 0.2, are graphed as functions of 9 in figure 7. The relative wavenumber Ki v for the 
nondissipative FCM (i.e., 5 = 0) is real-valued for 9 < 52° but has a nonzero imaginary part for 
9 > 52°. The relative wavenumbers Ki v for 5 = 0.1 and 0.2 have nonzero imaginary parts for all 
values of 9. Note that the real part of Ki V falls to zero at 6 = tt/2 in the absence of dissipation (i.e., 
5 = 0). 

Plots of the quantity w = Wi V , corresponding to the relative wavenumber Ki V of figure 7, are 
provided in figure 8. The negative phase-velocity condition Wi v < is satisfied 

(i) for 9 < 52° when 5 = 0, 

(ii) for 9 < 76° when 5 = 0.1, and 
(Hi) for 9 < 38° when 5 = 0.2. 

4 Discussion and Conclusion 

In isotropic dielectric-magnetic mediums, plane waves can propagate with phase velocity directed 
opposite to the direction of power flow under certain, rather restrictive, conditions [9]. However, 
the constitutive parameter space associated with anisotropic and bianisotropic mediums provides 
a wealth of opportunities for observing and exploiting negative phase-velocity behavior. Gen- 
eral conditions are established here for the phase velocity to be directed opposite to power flow 
for a particular class of bianisotropic mediums, namely Faraday chiral mediums. The theory has 
been explored by means of a representative example of FCMs, arising from the homogenization 
of an isotropic chiral medium and a magnetically biased ferrite. For our representative example, 
the negative phase-velocity conditions have been found to hold for propagation in arbitrary direc- 
tions — for both nondissipative and dissipative FCMs — provided that the gyrotropic parameter of 
the ferrite component medium is sufficiently large compared with the corresponding nongyrotropic 
permeability parameters. 

Previous studies [1]-[13] have emphasized the importance of the signs of constitutive (scalar) 
parameters in establishing the conditions for negative phase-velocity propagation in homogeneous 
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mediums. 2 In the absence of dissipation, negative phase-velocity propagation has been predicted 
in 

(i) isotropic dielectric-magnetic mediums provided that both the permittivity and permeability 
scalars are negative [3], and 

(ii) uniaxial dielectric-magnetic mediums when only one of the four constitutive scalars is negative 
[10]. 

Also, the conditions for negative phase-velocity propagation may be fulfilled by dissipative isotropic 
dielectric-magnetic mediums when only one of the two constitutive scalars has a negative real part 
[9]. The present study demonstrates that the condition for negative phase-velocity propagation can 
be satisfied by nondissipative FCMs with constitutive scalars that are all positive. Furthermore, 
these conditions continue to be satisfied after the introduction a small amount of dissipation. 

For the particular case of propagation parallel to the ferrite biasing field, the components of 
the time-averaged Poynting vector are null-valued in directions perpendicular to the propagation 
direction. In contrast, for general propagation directions, the time-averaged Poynting vector has 
nonzero components perpendicular to the direction of propagation. Further studies are required to 
explore the consequences of the negative phase-velocity condition ftR u • P(r) < for such general 
propagation directions. 

To conclude, more general bianisotropic mediums, particularly those developed as HCMs based 
on nonspherical particulate components, offer exciting prospects for future studies of negative 
phase-velocity propagation. 
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Figure 2: Bruggeman estimates of e HCM , £ HCM and fj, HCM as functions of The constitutive 
parameters of the component mediums are the same as in figure 1, but with /A 6 [0, 4] and / a = 0.35. 
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Figure 3: Calculated values of the relative wavenumbers Ki-i V as functions of /„, when 9 = and 
ix h g = 2, 3 and 4. The constitutive parameters of the component mediums are: e a = 3.2, £ a = 2.4, 
/i a = 2, e b = 2.2, = 3.5, and /j, b z = 1. Key: values corresponding to ^ = 2,3 and 4 

are represented by the solid lines, dashed lines, and broken dashed lines, respectively. Heavy lines 
indicate those relative wavenumbers which have nonzero imaginary parts; the real parts of such 
complex-valued relative wavenumbers are plotted. 
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Figure 4: Calculated values of Wi-i v as functions of / a , when 9 = and //* = 2,3 and 4. The 
constitutive parameters of the component mediums are the same as in figure 3. Key: Wi-i v values 
corresponding to fi b g = 2, 3 and 4 are represented by the solid lines, dashed lines, and broken dashed 
lines, respectively. Heavy lines indicate those w values that devolve from relative wavenumbers k 
with nonzero imaginary parts. 
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Figure 7: Real parts of relative wavenumbers Ki v as functions of 9 (in degrees) for the dissipation 
parameter b = 0,0.1 and 0.2 when f a = 0.35. The constitutive parameters of the component 
mediums are: e a = 3.2, = 2.4, /x a = 2, e b = 2.2 + ib, n b = 3.5 + i b, n\ = 1 + i 0.56, and 
fig = 4 + ilb. Key: Kj„ values corresponding to (5 = 0,0.1 and 0.2 are represented by the solid 
lines, dashed lines, and broken dashed lines, respectively. The heavy line on the graph for 6 = 
indicates those relative wavenumbers Ki v which have nonzero imaginary parts. 
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Figure 8: Calculated values of Wi v as functions of 6 (in degrees) for the dissipation parameter 
b = 0, 0.1 and 0.2 when f a = 0.35. The constitutive parameters of the component mediums are as 
in figure 7. Key: Wi V values corresponding to b = 0,0.1 and 0.2 are represented by the solid lines, 
dashed lines, and broken dashed lines, respectively. The heavy line on the graph for b = indicates 
those Wi v values which devolve from relative wavenumbers Kj„ with nonzero imaginary parts. 
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